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Abstract 

We present a systematic analysis of the angular distribution of -B — )• K*{-^ K'jt)1~^1~ decays with I = e,ii 
in the low recoil region (i.e. at high dilepton invariant masses of the order of the mass of the 6-quark) 
to account model-independently for CP violation beyond the Standard Model, working to next-to-leading 
order QCD. Prom the employed heavy quark effective theory framework we identify the key CP observables 
with reduced hadronic uncertainties. Since some of the CP asymmetries are CP-odd they can be measured 
without B-flavour tagging. This is particularly beneficial for Bs,Bs — )• <^(^ K~^K~)l~^l~ decays, which are 
not self-tagging, and we work out the corresponding time-integrated CP asymmetries. Presently available 
experimental constraints allow the proposed CP asymmetries to be sizeable, up to values of the order 0.2, 
while the corresponding Standard Model values receive a strong parametric suppression at the level of 
0(10"^). Furthermore, we work out the allowed ranges of the short-distance (Wilson) coefficients Cg in 
the presence of CP violation beyond the Standard Model but no further Dirac structures. We find the 
Bg — )• branching ratio to be below 9 x 10~® (at 95% CL). Possibilities to check the performance of 

the theoretical low recoil framework are pointed out. 



I. INTRODUCTION 



The exclusive rare flavour changing neutral current (FCNC) decay B — )• K*{-^ Ktt)1~^1~ with 
I = e, fi has high sensitivity to physics beyond the Standard Model (BSM) due to the large number 
of complementary measurements possible from the full angular distribution [Ij. Many works have 
focussed on the region of low dilepton invariant mass squared, q'^, typically taken within the range 
1-6 GeV'^. The latter is accessible to QCD factorisation OS], which has enabled systematic studies 
of CP-averaged observables as well as CP-asymmetries [IHH] . Intermediate values of fall into the 
narrow-resonance region dominated by the pronounced cc-resonance background from the decays 
B —7- K* {J/ip,tlj'} —7- K*l^l~ , recently studied in |9j including also low tails. At larger dilepton 
masses, at about q^ > 14 GeV^, follows the broad resonance region. The latter is characterised 
by the low recoil of the hadronic system. Here, the large values of ~ m^, where nib denotes 
the mass of the 6-quark, allow to perform an operator product expansion (OPE) |10l [TT] which, 
when combined with heavy quark effective theory (HQET) and the corresponding heavy quark 
form factor relations [12J, leads to powerful predictions, see [TO] and Hurth and Wyler in |13| . 
In fact, it has been shown recently that the heavy quark framework applied to the low recoil region 
results in a very simple amplitude structure of the decays B — t- K*{^ K'it)1^1^ [14J. Specifically, 
in the heavy quark limit, all three participating transversity amplitudes obey 

oc C^'^ X , i=^,||,0, (1.1) 

hence factorise into universal short-distance coefficients C^'^ and form factor coefficients /j. This 
feature can be greatly exploited to enhance the BSM sensitivity, to test form factor predictions 
against data and to check the goodness of the OPE framework. More explicit, the angular distri- 
bution of — 7- K*{^ Ktt)1^1^ decays allows for observables with the following salient properties, 
see [2] for details: 

i) The observable H!jl^ = 1 does not depend on short-distance coefficients nor on form factors. 

(2 3) (2) (3) 

ii) The observables Hj. ' depend on the short-distance coefficients only, and obey Hf, = Hj. . 
in) Several observables can be formed which depend on the form factors only. 

iv) The angular observables J7,8,9 5 which are odd under naive time-reversal, vanish. 

Beyond zeroth order in l/nih, the influence of the power corrections is weak because the Aqcd/?^6 
corrections are parametrically suppressed: The ones to the form factor relations enter with a 
suppression by small ratios of Wilson coefficients and the ones from subleading operators in the 
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OPE arise at 0{as) only. Moreover, the relevant hadronic matrix elements from both sources 
are not independent jlOI. While the latter matrix elements are currently not known from first 
principles for B — )• K*, model estimates suggest that they are at least not enhanced beyond the 
naive expectations |12j . 

In this paper we extend previous works [14J on the low recoil region by allowing for BSM CP 
violation. We work to next-to-leading order (NLO) in QCD and to lowest order (LO) in l/m^. The 
©(AqcD/'TT'b) corrections are taken into account in the estimation of the uncertainties. Further 
higher order corrections, including charm loops with gluons are power-suppressed at low recoil 
[TT] and not considered given the targeted precision. The consistency between the outcome 
of an analysis excluding and using only the high-g^ region data [H] supports the employed OPE 
framework. 

We propose and study CP observables with only subleading form factor uncertainties in Sections 



II A and II B , In Section II C we calculate mixing- induced time-integrated CP asymmetries relevant 



for the decays Bg, Bg — s- </>(— K^K ) l^l . In Section II D we give the relations between the CP 



observables and the angular distributions in i? — )• K*l^l or likewise Bg — )• (pl'^l decays. We work 



out the constraints on the complex-valued short-distance coefficients in Section III and summarize 
in Section |IV[ In an appendix we present the method used to estimate the uncertainties from the 
©(AqcD/'n-b) corrections. 

II. LOW RECOIL CP ASYMMETRIES 

Our aim is to extend our previous study of IS.B = 1 radiative and semileptonic decays [13] in the 
presence of CP violation. We use a model-independent framework with an effective Hamiltonian 

+ Cjmb [sa'"'PRb]F^,^ + h.c. + . . . , (2.1) 



where the ellipses denote contributions which we assume to be SM-like because they are either 
subdominant in the radiative and semileptonic 6 — )• s decay amplitudes or induced in the SM at 
tree level. With CP violation beyond the SM, the Wilson coefficients C7, Cg and Cio are complex- 
valued. All other Wilson coefficients are assumed to be SM-like, and are real-valued after factoring 



out the Cabibbo-Kobayashi-Maskawa (CKM) factors Vik, similar to Eq. (2.1). For details and 



definitions we refer to [TSj, which we follow closely. In particular, we take all numerical input as in 

h0.013 
-0.027' 



|14j except for the CKM one, which we calculate from the Wolfenstein parameters A = 0.812]'''^''^^^ 



A = 0.22543 ± 0.00077, p = 0.144 ± 0.025 and f) = 0.3421°;° [T5]. In the following we understand 
all Wilson coefficients to be evaluated at the scale ~ rrn,. The SM values of the most important 
ones are approximately, to next-to-next-to-leading logarithmic (NNLL) order, 

Cf^^ = _0.3, Cf^ = 4.2, Cfo^ = -4.2. (2.2) 



A. Decay amplitudes with CP violation at low recoil 



In a previous work [ll] we identified 

Pi = \{\C^? + \C'?) 



P2 



Cf + 



+ IC 



•R\2 



C7^|2) =Re<| (cf + J^Cf 



lol ' 



a 



10 f ' 



(2.3) 
(2.4) 



as the only independent short-distance factors which enter the observables of the angular distribu- 
tion of — )• K*{-^ Kit)1^1^ decays at low hadronic recoil. Here, s = q'^/m'^ ^-nd rhb = mh/niB, 
where denotes the mass of the B meson. The factor k accounts for the relation between 
the dipole (Ti^2,3) and (axial-) vector (y, j4i^2) form factors that can be calculated systematically 
|10t I12j. At lowest order in 1/mb and including 0{as) corrections, it reads 



The effective coefficients are written as 



Ceff p 



C3 + ^C4 + 20C5 + yC6 



+ 



OLs_ 



{C,-6C,)A{q')-C,FP{q'' 



and 



(2.5) 



(2.6) 



Cf = Cg + h{0, q^ 
1 



4^ ^ 11^ 2 32 

- + C2 + — C3 — - C4 + 52 C5 1^ Cg 



+ 



+ 



7C3 + -C4 + 76C5 + 



64 



+ 



Ci (i3((?2) + 4C(g2)) - 3C2 {2B{q') - C{q^)) - C^F^ (q^^ 



m 



7.^1 + 1^2] (1 + A„) + 2C3 + 2OQ 



(2.7) 



where we extended previous works |10| [T3] by including the doubly Cabibbo-suppressed contribu- 
tion proportional to = ^feKTs/C^b^s)- The latter is responsible for CP violation in the SM and 
appears only in the coefficient Cg^ with rn^/q^ suppression. We refer to |10|, I14j for more details 
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concerning the (real-valued) Wilson coefficients Cj<g as well as the LO and NLO QCD corrections 

encoded in the functions h(mi,q'^) and A,B,C, Fg'^'^\ respectively. 

In the presence of CP violation, there are four independent short-distance factors 

pi,P2 and pi,p2, (2.8) 

where the barred factors are obtained from the unbarred ones by complex conjugation of the Wilson 
coefficients Cj and the CKM factor A^. The building blocks describing CP violation are hence 

^Pi = Pi — Pi with i = l,2. (2-9) 



They can be written as 



using 



=4Imy-Im|c9 + K^C7 + y9^"U„| , (2.10) 
Ap2 = 2Imy-ImCio> (2-11) 



Y = Y^ + J^Yj (2.12) 



and the decomposition of Eqs. (2.6) and (2.7) into 



Cf = Ct + Yj, Cf = Cg + yg + KYt^ . (2.13) 

Note that Yg^"^ is real-valued, and we suppress the g^-dependence in the effective coefficients and 



the Yi throughout this work. It follows from Eqs. (2.10)-(2.11) that Api probes the weak phases 
of C7 and Cg, whereas Ap2 probes the weak phase of C^g. The imaginary parts of the Yi give rise 
to the strong phases and hence drive the magnitude of CP violation. 

In Fig. [T]we show the imaginary part of Yg (dashed curve) and Y (solid curve) from the OPE |10j 



at NLO QCD using Eqs. (2.6) and (2.7). The NLO QCD corrections to both luiYj^g are sizeable 
and lead to a reduction of the strong phases compared to the LO value of Im Y (dashed-dotted 
curve). Since Iml7 vanishes at LO, the NLO corrections constitute the leading contribution to 
this quantity. 

Also shown in Fig. [T]is the absorptive part Iml"'^'^ (dotted curve) obtained from a phenomenological 
fit to e~^e^ —7- hadrons data assuming factorization [16]. The Breit-Wigner amplitude matches the 
charmonium peaks of the branching ratios B{B — )• K*cc) for cc = J/^' and ^' . Note that for NNLL 
values of the Wilson coefficients C-^ and C2 and with present day data [T7| no fudge factor [TH] is 
needed. We apply this ansatz for the higher cc-resonances as well, where presently no i?-decay data 
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FIG. 1: The imaginary parts of Y (solid, blue) and Yg (dashed, red) in the OPE [TU] including NLO a^- 
corrections as functions of in the low recoil region. The LO result, where Imy = ImYg is shown by the 
dashed-dotted (orange) curve. The cc-resonance curve (dotted, black) shows the imaginary part of Yg from 
e~^e~' — !■ hadrons data |16l 117) . 

exist. The fit exhibits the local charm resonance structure from B — )• K*{cc) — )• K*l~^l~ decays. 
As can be seen in Fig. [T| the resonance contribution is of the same order as the OPE prediction at 
LO QCD and indicates comparable results after integrating over a sufficiently large region in the 
dilepton mass. However, we find a factor of ^■S between both approaches when using NLO QCD 
corrections and integrating over the low recoil region. 



B. The low recoil CP asymmetries 



At low recoil and within our framework (LO in l/mt, SM operator basis Eq. (2.1 ) ) all i? — )• K*l^l~ 
observables are, as far as short-distance physics is concerned, either proportional to pi, P2I P\ or 
short-distance insensitive [14j. Consequently, there are only two types of CP asymmetries: 

P2 _ P2 

CP ' CP \ } 



It is advantageous to define further 



rn^^Pi^ (2.15) 



('12") 

which is not independent of a^p . Since {p\ + pi) is positive definite, see Eq. (2.3), it is in general 



(2) 

better suited for normalisation then the denominator of a^p, which might cross zero and could 
make the theoretical uncertainties blow up. We note that a[|p equals the direct CP asymmetry in 



the rate, ^cp 



r-r 
r+r" 



Furthermore, a[?p equals Ap^ = ^b+Apb jjgj^ ^j^g Qp asymmetry of the 



(3) 

forward-backward asymmetry, whereas a^p corresponds to the low recoil transversity observables 



(2,3) 



introduced in Ref. 



No further CP asymmetries can be formed from the decays at low 



recoil beyond Eqs. (2.14)-(2.15) unless one considers neutral meson mixing, which we do in the 



(1 2) 

next section. Note that a^p are related to observables which require i3-flavour tagging, whereas 
can be extracted from untagged B meson samples [cf. Section 



II D 



Using A Pi <^ Pi, Pi the expressions for the CP asymmetries simplify to 



(1) _ 



(2) _ Ap2 



''CP 



2pi ' 2p2 - - pi 

The SM values of the CP asymmetries are induced at the order {m^/m^) ImA^ ~ 10^^ and are tiny. 
Since the CP asymmetries at low recoil are T-even only, a finite strong phase is needed for a finite 
CP asymmetry. The strong phase is roughly given as arg(y), yielding an additional suppression 
by another order of magnitude. Therefore, at low recoil 

,(1,2,3) 



^CP ' 



(3) _ Ap2 



^CP 



(2.16) 



^CP 



SM 



< 10- 



(2.17) 



Given the foreseen experimental precision, the CP asymmetries in the SM are therefore completely 
negligible due to their strong parametric suppression. Hence, any observed finite CP asymmetry 
is a signal of physics beyond the SM. 

Beyond the SM, the CP asymmetries at low recoil can be significantly enhanced. To estimate the 
order of magnitude of a^p ' we assume that jCy^g^iol are close to their respective SM values where 
ICy^l ^ ICg^l. Then, using Eq. (2.16) and ImY as shown in Fig. [l| we obtain, roughly. 



(1) 2 Imy 

"CP - If 12 I \n 12^ 



Im ( Cq + K 



O(O.l) 



2) Imyim C 



"CP ~ 



(3) 
"CP — 



10 



Re CgCjfo 



ag) < O(O.l) , 



Imyim C 



^<o(o.i) 



(2.18) 
(2.19) 
(2.20) 
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in agreement with |19H21j . Note that for very small BSM values of ReCgCJ'g the values of a^p 
become unconstrained. Furthermore, for ImC^g = the following relations hold 



*CP 



^CP ' 



^CP 



0. 



(2.21) 



To investigate more quantitatively the above CP asymmetries we define a BSM benchmark point 



-0.3i . 



Cg = +4.2 i, 



■-10 



-4.2 i, 



(2.22) 
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LO 


NLO 


cc-resonances 




+0.208 


+0.036 
-0.050 


SD 


±0.006 


SL 


+0.074 


+0.029 
-0.037 


SD 


±0.011 


SL 


+0.17 




+0.071 


+0.016 
-0.015 


SD 


±0.009 


SL 


+0.020 


+0.005 
-0.008 


SD 


±0.006 


SL 


+0.06 




-0.257 


+0.056 
-0.035 


SD 


±0.009 


SL 


-0.090 


+0.044 
-0.031 


SD 


±0.010 


SL 


-0.20 



TABLE I: The -integrated low recoil CP asymmetries at LO and NLO QCD for the BSM benchmark 



point (cf. Eq. (2.22|) with the main uncertainties from the variation of the renormahsation scale /if, (SD) 
and subleading power corrections (SL). Also given are the asymmetries using cc- resonance data [161 117) . 



which passes all the current experimental constraints. In particular, both interference terms 
ReCgQ andReCgC^ which are probed by ^FB(i? ^ {Xs,K*}l^r) and the ^ {X,, } /+/" 
branching ratios, respectively, are SM-like. Due to the maximal phases the benchmark values in- 
duce large BSM CP violation. The CP asymmetries evaluated at the benchmark point are given 
in Table |lj Throughout this work we use (..) to denote the integrated observables formed out of 
integrated angular coefficients following Ref. For the low recoil integration region we take 



14GeV^ < < 19.2 GeV^. We find that the main parametric uncertainty in (a[^p) stems from 
the variation of the renormahsation scale G [2.1, 8.4] GeV (SD). When varying the scale ^f,, the 



model-independent BSM contributions to the Wilson coefficients Eq. ( 2.22 ) are assumed to be given 
at the reference scale /x = 4.2 GeV and their renormahsation group evolution to /if, is taken into 
account in leading logarithmic approximation. While an efficient cancellation of the /ife-dependence 
is at work in the sum Cf^ = Ci + Yi, the numerators of the CP asymmetries depend on the product 
of ImCj and ImY, and result in the large reported uncertainties. The subleading corrections to the 
Isgur-Wise form factor relations and transversity amplitudes (SL) constitute another major source 
of uncertainty. Its estimate is explained in Appendix [A} 

As can be seen from Table |I] the impact of the NLO corrections is sizable on the CP asymmetries. 
The LO predictions are about a factor 3 larger than the NLO ones due to large destructive NLO 
contributions to Iml". In fact, concerning Iml^ the NLO corrections constitute the leading contri- 
bution which also implies a large scale uncertainty at NLO, but the NLO corrections are sizeable 
in Imlg, too. 

Also shown in Table |I] are the CP asymmetries calculated using a phenomenological ansatz with 
cc-resonances |16t [T7] . They are in the general ballpark of the OPE ones, between the LO and 
NLO findings, and somewhat smaller than the LO results. 
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C. Untagged CP asymmetries with meson mixing 



We consider the decays Bg^Bg — )• K~)l~^l~ which especially for muons are of great im- 

portance for hadron collider experiments. We follow closely [3] to which we refer for details on the 
full angular distribution [1]. 

To account for neutral meson mixing, time-dependent transversity amplitudes need to be intro- 
duced: 

ALJR{t) ^ A'^/^'iBsit) ^ K+K-)J+l-) , 

A^/'^it) = ^ K+K-)J+l-) , (2.23) 

T / f? — T / f? — 

where Aa (t), {Aa (t)) denotes the amplitude for a meson born at time t = as a Bg, {Bg) 
decaying through the transversity amplitude a =_L, ||, at later times t. 

For the time evolution the following parameters which involve the un-mixed amplitudes at t = 
play an important role 

^L/R ^ -^"^""(Q) , (2.24) 

A^^^'mSw ^ -dw] 



where [6w — )• —6w] implies the conjugation of all weak phases in the denominator. Here denotes 

>SM = 2arg(F,*1 



the phase of the Bg — Bg mixing amplitude which is very small in the SM, = 2 arg(V^* VJ;,). 



The untagged rates dT + dV can then be written as |22j 

A^{t)Alit) + A^{t)Alit) = IaMAUO) (2.25) 

X [{l + Var]bCaCb){e-^'^' + e-^"') + {r]aU + rjhCb){e-^''-e-^"')] , 

where the chirality indices L,R are suppressed for brevity. Here, r;o,|| = +1 V± — ~1 
the CP eigenvalues of the final state and ^l{H) denotes the width of the lighter (heavier) mass 
eigenstate of the Bg system. We also neglect CP violation in mixing, which is bounded by the 
semileptonic asymmetry for Bg mesons I^slI ^ 0(10"^) |23j . 
After time- integration follows from Eq. ( 2.25| ) 



/ dt [A^miit) + Aa{t)Al{t)] = ; / [1 + VaVbUl - y iVa^a + VbCb)] , (2-26) 

Jo - y ) 

where P = (Pl + P_f/)/2 and AP = P^ — Th denote the average width and the width difference, 
respectively, and y = AP/(2P). 
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Due to the simple transversity structure at low recoil [13] , see Eq. (1.1), there are only two different 



time evolution parameters, and ^r, universal for all the _L, ||,0 amplitudes. We obtain 



e^/^ = e-'*M ^oTCio + ^'T" C. + Y + KY^'"'^ ^^^^^^ 

CI T Ci*o + q + y + xiYp 



In the absence of strong phases we find |Cl/_r| = 1 and in the absence of CP violation in the 
rare decays holds = 1 as well. In the SM, CP violation in these parameters is very small: 

ieL//?|-l = 0(K/m2)ImA„). 

CP-odd observables allow to measure CP violation without tagging the flavour of the initial B 
meson (if the asymmetry between Bg and Bg production is known). Among the available coefficients 
in the angular distribution four of them, J5,6,8,9! are CP-odd [Illl]- However, in the low recoil region 
the coefficients ^7,8,9 vanish [Hj, leaving only J^^ for untagged CP measurements. 



For J5 we obtain at low recoil from Eq. (2.26): 



dme[A^{t)A]:*{t)+A^{t)A]:*{t) - (L ^ R)] = foh x f^J^ (2-28) 



with 



Amix = 2p2(iar + le^r - 2) + pide^r - lan ■ (2.29) 

The formula for Jg is identical after changing the transversity index to ||. 

In order to reduce non-perturbative uncertainties, we choose combinations of the following un- 
tagged, time-integrated quantities for normalisation: 

n. = Hdi [\AHt)\' + \AHt)\' + iL^ R)] = ff X ~ 2r;.^ C^..) ^ ^^ ^^^ 

Jo i u ~ y J 

which can be obtained from the angular observables ^1,2,3- Here, 

Brm. = PidaP + IChP + 2) + 2p2(|eRp - lap) , (2.31) 

Crmx = PiRe{^L + (r} + 2p2Re{^B, - ^l} • (2.32) 

Note that Cmix only depends on the mixing phase $Af. 



(3) 

Normalizing Eq. (2.28) to y^njn^ yields a mixing-induced analogue of a^p.: 



a^if = ^ . (2.33) 



Note that at low recoil a^p^ is insensitive to the sign of y. Simultaneously, the sensitivity to 
is very low since it enters via Cmix only. In the limit y — t- holds 

^mix ^ ^mix ^ _2 jPl + 4/92)P2 " 2piP2Pl + (Pl "4^1)^2 
\Bmix\ {pl+4.pl)Pi-SpiP2p2 + {p\-^pl)pi' 
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For A/9j <^ Pi, Pi this simplifies further to 

a^r = 4p- (2-35) 

The use of the coefficient Jg with normahzation to y/n^^nj^ leads to a second possibility to measure 
the very same asymmetry a^^. 

The sensitivity to the i3<j-mixing parameters y and is low for realistic values of y < 0(0. 1). 
We find for the g^-integrated asymmetries that, model- independently, |(a^i,^)/(a[?^) - 1| is below 
a few percent and hence unlikely to be measured in the foreseen future. 

D. CP asymmetries from the angular distribution 

Here we summarise the relations between the low recoil CP asymmetries a^p ' and the CP asym- 
metries studied previously in the literature. We begin by showing various possibilities to extract 
the ajjp from the angular distribution. We neglect the small corrections from finite lepton masses 



mi in kinematical factors /3; = — Amj/q^. 

As already mentioned, a[lp equals the total rate asymmetry ^cp given as 



r r 

= AcP = . (2.36) 



T + T 

Here, form factor uncertainties cancel at low recoil. 

(2) 

The asymmetry Oq^ can be extracted in a multitude of ways from the ratios 



«CP= ; , o Ra,= ^^, Ra, = RaM^^J^], (2-37) 



^(2) _ Rgb — Rgb ^ _ Z^g '^a^a 

Rab + Rab ^b'^bJb' 

which provides a cancellation of the form factor uncertainties in the low recoil region. The and 

"-" sign in the third relation apply for CP-even i = 1, 2, 3, 4 and CP-odd i = 5, 6 coefficients Jj 

in the angular decay distribution, respectively. The ratios Rab can have a = 5, 6 and 6 = 1, 2, 3, 4 

and the Cj are auxiliary real numbers. For example, Ap^ can be recovered for CaJa = Je and 

Y.b'^hJb = Ji- '/2/3 (with Ji_2 = 2Jf 2 + Jf,2) 

„(2) _ .CP _ ^FB + ^FB . 
«CP - ^FB - _ • {1.6^} 

The CP asymmeries Ai = 2{Ji - Ji)/d{T + f)/dq^ and A^ = -2{Ji - Ji)/d{T + f)/dq^ defined in 
[3] are related in the low recoil region to the CP asymmetries a^p as 

, _ (1) 3 fl - f\\ aD - _ (1) _SL /o/|| 

AD - _„(3) V A /o/-^ A _ A3) o /-L/jl oQ^ 

A - «CPX ^^2 + ,2 + ^2' ^6-acpX3^2 + ,2 + ^2 ' (2.39) 
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whereas Aj^^g vanish. Unhke in a^p^'^'' in the asymmetries Eq. (2.39) the form factors do not drop 
out and the corresponding uncertainties do not cancel. 

(3) (2) (3) 

The CP asymmetry a^p is identical to the low recoil transversity observables Hj^ and Hj, 
when measuring them in an untagged sample containing an equal number of B and B mesons 



^CP 



J5 — J5 



-2(J| + J2^)[2(J| + J|) + (J3 + J3)] 

Je — Jq 



for K 



(2) 
T 



for H, 



(3) 



(2.40) 



^2^J4{Ji + J|Y-{Js + J3Y 
The mixing-induced, time- integrated CP asymmetries A^"^^^ and A™*^ defined in [2] are related 



to as follows 



A Dmix 
^5 



mix 
-Qqp X 



non± 



Al 



aJSh^ X 3 



n_Lnj| 



(2.41) 



"11 



where the rii have been introduced in Eq. (2.30). As in Eq. (2.39), and unlike in a^p^, the asym- 



metries given in Eq. (2.41) exhibit a residual form factor dependence. 

Ways to extract the angular coefficients Jj from the full differential decay distribution have been 
given in Ref. [1], to which we refer for details and the definitions of the kinematic angles 61^ Ok* 
and 4>. For example, requires two bins, cos6'; G [— l^O] and [0, 1], while J5 can be extracted from 
the single-differential distribution in the angle with additional bins, cos^x* £ ["1)0] s-nd [0, 1]. 



III. MODEL-INDEPENDENT AS = 1 CONSTRAINTS 



We perform a global analysis of the available h — )• sl^l^ decay data in the presence of BSM 
CP violation through the Wilson coefficients Cjqiq, i.e., allowing them to be complex- valued. 
We follow the same approach and the same data sources as presented in |14j to perform a six- 
dimensional scan of the magnitudes jC^g and phases </)7,9,io = argC^g iq. We use the following 
ranges and binning 

\Cj\ £ [0.30,0.35] AIC7I =0.01, (pr £ [0,2ir) A(l)j 

\Cg\ G [0, 15] A|Cg| = 0.25 , G [0, 2tt) A^g 

ICiol G [0, 15] A|Ciol = 0.25 , cj^io G [0, 27r) A(/>io 

The narrow range for jC^I ~ \Cy^\ is justified by the good agreement of the measured B — )• Xgj 
branching ratio with its SM prediction |24t I25j . For the scan we used and developed further 
EOS [26j, a program for the evaluation of flavour observables. 
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FIG. 2: The constraints on |Cg| and \Ciq\ from the experimental data as collected in [14^. The areas 
correspond to 68% CL (red) and 95% CL intervals (blue). In the left plot (a), data from the low recoil 
region has been excluded, while the right plot (b) has been obtained using the full set of available data. The 
(light green) square marks the SM. 

In order to visualize the constraints, we project the 68% and 95% confidence regions of the six- 
dimensional scan onto the |Cg|~|C]^ol plane, shown in Fig. [2] The projections onto |Cg|-(/)9, |C]^q|-(/>io 
and (j)g~(pio are obtained likewise, and are shown in Fig. [sj The data from B — ?• K*l'^l~ decays 
in the low recoil region provide powerful additional constraints as can be seen by comparing the 
results with or without including them. We find good agreement between the SM and the data. 
The scan procedure returns the allowed ranges, see Fig. [2j 

1.0<|C9|<6.5 (0.3 <|C9| < 7.3) , 

2.0 <|Ciol < 5.3 (1.5 <|Cio| < 5.8) , (3.2) 

at 68% CL (95% CL). Due to its smallness the above finite lower 95% CL-bound on |Cg| is sensitive 
to the discretisation of the scan, A|Cg| = 0.25, and is subject to corresponding uncertainties. 



From Eqs. (2.2) and (3.2) we find for the branching ratio of the decay Bg — s- with respect to 

its SM value a maximal enhancement by a factor 1.9. Employing for the decay constant of the Bs 
meson fss = 231(15)(4) MeV [27J, we obtain the 95% CL upper limit B{Bs ^ /U+zU") < 8 x 10"^. 
The corresponding SM value is given as B{Bs — s- /t^^^~)sM = (3.1 ± 0.6) x 10^^ with the dominant 
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Large Recoil + Inclusive Large Recoil + Inclusive Large Recoil + Inclusive 




FIG. 3: The constraints on complex- valued Cg and Ciq from the experimental data as collected in TT. The 
areas correspond to 68% CL (red) and 95% CL intervals (blue). The upper row shows the constraints using 
the large recoil region of the exclusive and inclusive decays only. In the lower row, all data from both the low 
and the large recoil regions are used. The solid (black) lines mark the 68% CL regions obtained by employing 
in addition a hypothetical measurement of (a|jp) — 0.074 ± 0.01 with a central value corresponding to the 
one at the BSM benchmark point (cf. Eq. ( 2.22 1). The (light green) square denotes the SM value of {Cq,Ciq), 



while the benchmark point is denoted by the black circle. 

uncertainty stemming from the decay constant. Using fss = 256(6) (6) MeV [28], we obtain a 
slightly larger upper limit B{Bs H^H') < 9 x lO^^, and B{Bs Ai+^~)sM = (3.8 ±0.4) x IQ-^. 
The SM region of Bg — )• fj.~^fi~ decays will be at least partially accessed by the LHCb experiment 
with the 2011-2012 LHC run with projected luminosity up to around 2 fb^^ (29] . 
We find the approximate 95% CL ranges for the phases, see Fig. [sj 

^<arg(C9Ci*o)< ^, (3.3) 

with corresponding 27r-periodic branches. We note that the experimental information entering our 
scans stems from CP-conserving data only and the constraints on the BSM phases are currently 
weak. While the S-factories already measured the rate asymmetry of — )• K*l^l~ decays at the 
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level of 0(0. 1) |23| . these constraints are not included in our analysis because they are given for 
the total integrated rate only, or are inappropriately binned such as the high rate asymmetry 
measurement by BaBar, ^cp = 0-09 it 0.21 it 0.02, given for > 10.24 GeV^ and excluding the 
^''-peak [30] . We checked explicitly that an (a^p) measurement at the level of the latter with 
theoretical uncertainties taken into account is not significant in the scan. We find the values of the 
CP-observables (a^p ) to be within -0.2 and +0.2, while (o^p) unconstrained by current data. 
To illustrate the impact of a future measurement of the CP asymmetry, we add hypothetical data 
with a reduced experimental uncertainty (app) = 0.074 it 0.01 to the scan. The central value is 
inspired by the BSM benchmark, see Table |I| As can be seen from Fig. [3| the extended data set 
(solid, black lines) adds complementary constraints on the phases, which become challenging to 
the SM (green square). 



IV. SUMMARY 

Building on previous works on CP symmetries [H] we identified CP asymmetries a^p, i = 1,2,3 
and with no leading form factor dependence from the angular distribution of i? — )• X*(— >■ 

Ktt)1^1~ and of Bs,Bs — t- <^(— t- K^K^)l^l^ decays at low hadronic recoil. The simple amplitude 
structure following from the heavy quark framework of Ref. jlO] in this kinematical region has been 
crucial in doing so. We find that the largest uncertainty in these CP asymmetries stems from the 
renormalisation scale dependence at NLO in as, which is sizeable, followed by subleading 1/mf, 
corrections, see Section [TlBj 

Being strongly parametrically suppressed in the SM, the CP asymmetries are nulltests of the 
SM. At the same time the available experimental constraints allow for large BSM effects in the 
(low recoil (/^-integrated) asymmetries (ocp^^) up to ~ 0.2. The asymmetry (a^p) is due to its 
possibly vanishing normalisation presently unconstrained. The mixing asymmetry in Bs,Bs — >• 
0(— 7- K^K^)l^l~ , {a^Qp) exhibits for realistic values of the mixing parameters little sensitivity to 
the latter and is numerically close to (a^p)- Note that both o^p^ and a^p do not require flavour 
tagging. 



By testing the effective theory Eq. (2.1 ) against the existing data of 6 — )• 57 and h — )• sl^l~ decays 
in the presence of BSM CP violation, we extract the allowed ranges of the Wilson coefficients 
Cg shown in Figs. [2] and [3| We find consistency with a related recent analysis for real- valued 
coefficients [14J and with the SM. Parameter points with order one deviations from the SM are 
presently allowed. 
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To maximize the exploitation of data we strongly suggest to provide the future CP symmetries 
and asymmetries in g^-bins accessible to systematic theory calculations, such as 1-6 GeV^ and 
> 14GeV^, similar to the common binning used in both recent Belle and CDF analyses 132] . 
For completeness we give here the low recoil SM predictions for the basic CP-averaged observables, 
the branching ratio, the forward-backward asymmetry and the fraction of longitudinally polarised 
K* mesons, 

/•(m_g-mx*)2 
' 14 GeV^ 



lO'^ X 



dq 



;d5sM 



dg2 

(^fb)sm 



2.86 



hO.87 
-0.74 



+0.11 




fO.lO 


+0.08 




FF -0-1° 


SL 


-0.19 


CKM 


SD 



-0.41 ±0.07 



±0.007 



FF 



SL 



+0.002 
-0.003 



SD 



SM 



0.35 



h0.04 
-0.05 



FF 



±0.003 



SL 



(4.1) 
(4.2) 
(4.3) 



respectively. The predictions are based on the improved uncertainty estimate for the subleading 
power corrections of Appendix [A| and updated CKM input from |15] . but follow |14j otherwise. 
The largest uncertainty in the above observables stems from the form factors (FF). Uncertainties 
smaller than a permille are not given explicitly. 

We stress that the transversity observables allow for consistency checks of the theoretical low recoil 
framework. Since the OPE-breaking corrections generically will spoil the transversity amplitude 



relation Eq. (1.1) on which the predictions i) - iv) listed in the Introduction are based on, the 



performance of the employed heavy quark framework can be tested experimentally. 
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Appendix A: Estimate of the subleading power corrections 

The subleading power corrections to i? — t- K*l~^l~ decays at low recoil arise from the form factor 
(Isgur-Wise) relations beyond the heavy quark limit and the contributions of subleading opera- 
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tors to the OPE [TU]. Both contributions involve the same (three) HQET form factors, which 
are essentiahy unknown (see [T2j for model estimates) but in principle are accessible to lattice 
calculations. 

The subleading corrections to the form factor relations enter the decay amplitudes multiplying the 
small coefficient Cj. The subleading OPE-corrections involve presently unknown Wilson coefficients 
of order a^. Their knowledge would require a generalisation of the 2-loop calculation of [33J for 
off-shell quark states. The latter gives in general complex-valued results, introducing new strong 
phases. Both corrections can therefore be parameterised as 

Aqcd , „ /,A 



fi ~ ±- 



mb 



i = a,b,c. 



(Al) 



The exact form of the fj can be inferred from fj = Cf, with the Tj given in Ref. |10j . As already 
noted, the strong phases 6i are currently not known. 

The transversity amplitudes depend on the fj as follows, see [13] for details. 



(Cf tCio) + '^— Cf + r„ 



A 



L,R 



(A2) 
(A3) 



A 



L,R 



(Cf =F Ciq) + K—z-Cj^ 



fo 



iNiriB 



(1 



mj^,){l + rhK^YhAi - XfcA 



2 rriK* (1 + mK-'^y s 



(A4) 



which mildly breaks the universality of Eq. (1.1). In the numerical implementation we vary the 
fi and 5i for i = a,b,c independently within |fj| < 0.1 and 6i £ [— 7r/2, -|-7r/2] and allow for both 



signs in Eq. (Al). The resulting uncertainty is termed (SL) in this work. The SM predictions of 



some basic CP-conserving observables are given in Section |IV[ Note that this estimation improves 
on our previous works ^4j, where we introduced for each of the transversity amplitudes one real 
scaling factor for the corrections to the form factor relations and additionally for each left- and 
right-handed amplitude six real scaling factors in order to estimate the subleading OPE-corrections. 



[1] F. Kruger, L. M. Sehgal, N. Sinha and R. Sinha, Phys. Rev. D 61, 114028 (2000) [Erratum-ibid. D 63, 

019901 (2001)], |arXiv:hep-ph/9907386| . 
[2] M. Beneke, T. Feldmann, D. Seidel, Nucl. Phys. B612, 25-58 (2001), |hep-ph/0106067| . 
[3] M. Beneke, T. .Feldmann, D. Seidel, Eur. Phys. J. C41, 173-188 (2005), |hep-ph/ 04 12400) . 
[4] C. Bobeth, G. HiUer and G. Piranishvili, JHEP 0807, 106 (2008), iarXiv:0805.2525l [hep-ph]]. 



17 



[5] U. Egede et at, JHEP 0811, 032 (2008), |arXiv:0807.2589| [hep-ph 



[6] W. Altmannshofer et al, JHEP 0901, 019 (2009), |arXiv:0811.1214| [hep-ph]] 
[7] C. Bobeth, G. Hiller and G. Piranishvili, |arXiv:0911.4054| [hep-ph]. 



[8] U. Egede, T. Hurth, J. Matias et al, JHEP 1010, 056 (2010), |arXiv: 1005.0571 ' [hep-ph]]. 
[9] A. Khodjamirian, T. Mannel, A.A. Pivovarov and Y.M. Wang, JHEP 1009 (2010) 089, jarXiv:1006.4945 
[hep-ph]]. 



[10] B. Grinstein and D. Pirjol, Phys. Rev. D 70, 114005 (2004), |arXiv:hep-ph/0404250 



[11] M. Beyhch, G. Buchalla and T. Feldmann, |arXiv:1101.5118| [hep-ph 



[12] B. Grinstein and D. Pirjol, Phys. Lett. B 533, 8 (2002), |arXiv:hep-ph/0201298] . 



[13] J. L. Hewett et al, |arXiv:hep-ph/0503261 



[14] C. Bobeth, G. Hiller and D. van Dyk, JHEP 1007, 098 (2010), [arXiv:1006.50T3 [hep-ph 



[15] J. Charles et al. [ CKMfitter Group Collaboration ] , Eur. Phys. J. C41, 1-131 (2005), | hep-ph/0406184| . 

We use the numerical results as presented at ICHEP2010, Paris, France, July 22-28, 2010. 

[16] F. Kruger and L. M. Sehgal, Phys. Lett. B 380, 199 (1996), |arXiv:hep-ph/9603237| . 

[17] K. Nakamura et al. [Particle Data Group], J. Phys. G 37, 075021 (2010). 

[18] A. Ah, P. BaU, L. T. Handoko and G. Hiller, Phys. Rev. D 61, 074024 (2000), |arXiv:hep-ph/9910221| . 

[19] G. Buchalla, G. Hiller and G. Isidori, Phys. Rev. D 63, 014015 (2000), |arXiv:hep-ph/0006136| . 

[20] F. Kruger and E. Lunghi, Phys. Rev. D 63, 014013 (2001), (arXiv:hep-ph/0008210| . 



[21] A. K. Alok, A. Datta, A. Dighe, M. Duraisamy, D. Ghosh and D. London, arXiv:1103.5344 [hep-ph]. 



[22] R. Fleischer and I. Dunietz, Phys. Rev. D 55, 259 (1997), |arXiv:hep-ph/9605220| . 



[23] E. Barberio et al. [Heavy Flavor Averaging Group], |arXiv:0808.1297| [hep-ex 



Online |http: / /www. slac. Stanford. edu/xorg/hfag| from March 2010. 
[24] M. Misiak et al, Phys. Rev. Lett. 98, 022002 (2007), |arXiv:hep-ph/0609232| . 
[25] T. Becher and M. Neubert, Phys. Rev. Lett. 98, 022003 (2007), |arXiv:hep-ph/0610067| . 



26] D. van Dyk et al, Online http://project.het.physik.tu-dortmund.de/source/eos/tag/lowrecoil-cp 



[27] E. Gamiz, C. T. H. Davies, G. P. Lepage, J. Shigemitsu and M. Wingate [HPQCD CoUaboration] , 



Phys. Rev. D 80 (2009) 014503, (arXiv:0902.1815| [hep-lat]] 



[28] J. Simone et al [Fermilab Lattice and MILC Cohaborations] , PoS LATTICE2010 (2010) 317. 
[29] M. Palutan on behalf of the LHCb Collaboration at BEAUTY 2011, 13th International Conference on 
B-Physics at Hadron Machines, April 4th-8th 2011, Amsterdam, The Netherlands. 



B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 102, 091803 (2009), |arXiv:0807.4119| [hep- 
ex]]. 



[31] J. T. Wei et al [BELLE Collaboration], Phys. Rev. Lett. 103, 171801 (2009), |arXiv:0904.0770 [hep-ex 



[32] T. Aahonen et al [CDF Collaboration], arXiv:1101.1028 [hep-ex]. 



[33] D. Seidel, Phys. Rev. D 70 (2004) 094038, |arXiv:hep-ph/0403185| . 



18 



